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Abstract 

In this paper we consider the Glauber dynamics for a disordered ferromagnetic Ising 

1 model, in the region of phase coexistence. It was conjectured several decades ago that the 
spin autocorrelation decays as a negative power of time |HF87j . We confirm this behavior 

p I by establishing a corresponding lower bound in any dimensions d J? 2, together with an 

^ upper bound when d — 2. Our approach is deeply connected to the Wulff construction for 

the dilute Ising model. We consider initial phase profiles with a reduced surface tension on 
their boundary and prove that, under mild conditions, those profiles are separated from the 
G (equilibrium) pure plus phase by an energy barrier. 
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1 Introduction and definitions 
1.1 Introduction 

For many years the Ising model and the corresponding Glauber dynamics have been a very 
active research field. In the 1990's the asymptotics of the spectral gap of the Ising model with 
free boundary condition were connected to surface tension, sec [Tho89, Mar99] and references 
therein. The inversion time of the infinite volume Ising model phase under a small field was 
then related to Wulff energies [Sch94. SS98]. In the last decade, precise estimates were achieved 
for the spectral gap and mixing time }BM02| . Recently impressive moves towards evidence of 
Lifshitz behavior and mean-curvature displacement of interfaces were achieved in [MT10, LS10, 
ILMST101 ICMST10] . 

The focus of the present paper is on the consequences of the presence of disorder on the 
dynamics, in the phase coexistence region. Since [Mar99j (and references therein) it has been 
known that dilution in the Ising model triggers slow, non-exponential relaxation to equilibrium 
in the Griffiths phase. Here we focus on the phase coexistence region, which means that, at 
equilibrium, the system can be either in the plus or the minus phase. This setting was considered 
already in [HF87J where heuristic discussions suggested that autocorrelation decays as a negative 
power of time. In the present paper we turn these heuristics into rigorous proofs. Previous stages 
of this project were the adaptation of the coarse graining and of the Wulff construction to the 
disordered setting, see |Wou081 rWoun9] respectively. 



Our main result is a lower bound on the autocorrelation (Theorem 2.2) which validates the 
heuristics of [HF87J. That is to say, when both the initial configuration includes a droplet of the 
minus phase and the surface tension on the boundary of the droplet is smaller than its quenched 
value, the system must cross an energy barrier before the droplet can disappear. Interestingly 
we show that the energy gap can be computed on continuous evolutions of the droplet, a slight 
improvement in comparison with the usual scheme of computing the bottleneck as the maximum 
gap in energy over all intermediate magnetization, as in [Mar99, BI04J. 

We also present in this work an upper bound on the autocorrelation when d = 2 (Theorem 



2.4) together with some consequences of our estimates on the typical spectral gap and mixing 



times in finite volume (Theorems 2.5) 



We would like to mention that, although is it the case here, we do not expect that dilution 
always slows down relaxation. Indeed, for the infinite volume dilute Ising system subject to a 
small positive external field, the dilution has a catalyst effect on the transition from the minus 
to the (equilibrium) plus phase, cf. [BGW12) . 

The organization of the paper is as follows: in the remaining part of the current Section, 
we define the dilute Ising model and the Glauber dynamics. We also introduce the necessary 
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tools and technical assumptions. Then in Section [2] we present our main results. Heuristics and 
proofs are given in Section [3} 



1.2 The dilute Ising model 

The canonical vectors of R, d are denoted by (ei)j=i...d. For any x = YH=i x i e i = • • • 
R d we consider the following norms: 




d 

\\x\\i = \xj\, \\x\\2 = ( ^~] x t 1 an d Hxlloo = max . (1-1) 
i=l \i=l / *~ 

Given x, y G Z d we say that x, y are nearest neighbors (which we denote x ~ y) if they are at 
Euclidean distance 1, i.e. if||x — y||2 = l- To any domain A C Z d we associate the edge sets 

E(A) = {{x,y} : x,y G A and x ~ y} (1.2) 
and E W (A) = {{x, y} : x G A, y € Z d and x ~ y} . (1.3) 

We consider in this paper the dilute Ising model on 7L d for d ^ 2. It is defined in two steps: first, 
the couplings between adjacent spins are represented by a random sequence J = {Je) e eE(w, d ) °f 
law P, such that the (Je) e eE(w, d ) are independent, identically distributed in [0, 1] under P. For 
convenience we write J = [0, l] s ( zd ) the set of possible realizations of J. Given Ac^a finite 
domain, J €L J and a spin configuration a G = {cr : Z d — )• {±1} : = 1, Vz ^ A}, we let 

#a» = - E ( L4 ) 

e={a;,i/}e£ !1 "(A) 

the Hamiltonian with plus boundary condition on A. The dilute Ising model on A with plus 
boundary condition, given a realization J of the couplings, is the probability measure fj,^ + on 
that satisfies 

Ma' + (M) = (-§ #a'») , Va G £+ (1.5) 



where /3 ^ is the inverse temperature and is the partition function 

P 



Consider 



m p = lim E//-+ (<r ) (1.7) 

the magnetization in the thermodynamic limit, where A^v is the symmetric box A^r = {— A r , . . . , A^} 
and IE the expectation associated with P. When mp > the boundary condition has an 
influence even if it is arbitrary far away from the origin. In particular the decreasing limit 
= limjv 4 /^' + ' does not coincide with = lirnzv t A 4 ^'"^ (P-almost surely). The two 

infinite volume measures E//~' + and E/io'~ are respectively the plus and minus phases. 



Z ti= E ex P (-|i^+(a)). (1.6) 
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It was shown in |ACCN87] that the dilute Ising model undergoes a phase transition at low 
temperature when the random interactions percolate. In our settings, this means that the critical 
inverse temperature 

/3 c = inf{/3^0:m^>0}, (1.8) 

which is never smaller than 0^ UTe - the critical inverse temperature for the pure Ising model 
corresponding to J = 1 - is finite if and only if P( J e > 0) > p c (d) where p c {d) is the threshold 
for bond percolation on 7L d . 

1.3 The Fortuin-Kasteleyn representation 

The Ising model has a percolation-like representation which is very convenient for formulating 
two of the fundamental concepts for the study of equilibrium phase coexistence: renormalization 
and surface tension. We call 

n = <jw : E(% d ) -> {0,1}} 

the set of cluster configurations on E{'E d ), and for any oj G Q and E C ECE d ) we call oj\ e the 
restriction of oj to E, defined by 

, , ( oj e if e G E 
^ E > e ~ { else. 

The set of cluster configurations on E is £l E = {u\ e ,oj G Given a parameter q ^ 1 and an 
inverse temperature j3 ^ 0, a realization of the random couplings J : E(T* d ) — > [0, 1], a finite 
edge set E C E(% d ) and a boundary condition 7r G fiec we consider the random cluster model 
^E^p 9 on defined by 

^;/(M) = -^IIPe e (l-Pe) 1 -^x^H, Va;eO £ (1.9) 

E,/3 e&E 

where p e = 1 — exp(— /3J e ), C e (oj) is the number of clusters of the set of vertices in 7L d attained 
by E under the wiring oj V ir such that (oj V 7r) e = max(w c , vr e ), and Z J E ,1T ^ q is the renormalization 
constant that makes of Q^a 1 a probability measure. 

For convenience we use the same notation for the probability measure ^ E n p 9 and for its 
expectation. Most of the time we will take either n = /, where / is the free boundary condition: 
f e = 0,Ve G E c , or tt = w where w is the wired boundary condition: w e = 1, Ve G E c . When 
the parameters q and /3 are clear from the context we omit them. Given 1Z a compact subset 
of R d (usually a rectangular parallelepiped) we denote by <&^ 7r the measure ^ E ^ n ^ on ^ ne 

cluster configurations on E(JZC\7L d ), where 1Z stands for the interior of 1Z. In particular, for any 
g, h : fl — > R the quantities ^^(g) and &^(h) are independent under P when 7£i n 72-2 = 0- 

The connection between the dilute Ising model and the random-cluster model was made 
explicit in [ES88J. Consider the joint probability measure 

= II (Pera-Pe) 1 ^, V(^)£E{x(l E(A) 

Z A,/3 ee£™(A) 

where p e = 1 — exp(— /3J e ), cr -< a; is the event that a and w are compatible, namely that 
oj e = 1 a x = a y ,ye = {x, y} G E W (A), and is the corresponding normalizing factor. 

Then, 
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i. The marginal of Vf^'lj on the variable a is the Ising model fi'^ + , 

ii. Its marginal on the variable u is the random-cluster model &"gfj^ p with wired boundary 
condition w and parameter q = 2. 

iii. Conditionally on uj, the spin <r of each connected component of A for uj (from now on cluster) 
is constant, and equal to +1 if the cluster is connected to A c . The spin of all clusters not 
touching A c are independent and equal to +1 with a probability 1/2. 

iv. Conditionally on a, the edges are open (i.e. uj e = 1 for e = {x, y}) independently, with 
respective probabilities Ve&a x ,ay 

Furthermore, the distribution ^^Q 9 increases with j3, J, tt, it satisfies the FKG inequality 
and the DLR equation, cf. [ACCN88 ] .' 

1.4 Slab percolation 

We say that slab percolation holds under E<3?^jf^ ^ when either 

d^3 and 3H £ N* inf inf E^ f ' g Jx y) > 0, 

or d = 2 and 3k : N* i— > N* lim E^c^' 9 «(3 a horizontal crossing for uj) > 

AM-oo 6 iV,K(iV)^ v to y 

where Sn,h = {1> • • • > -/V}^ -1 x {1, . . . , H} is the slab of height H. The critical threshold for 
slab percolation is 

/3 C = inf j/3 ^ : slab percolation occurs under E^jf^ ^ j . (1-10) 

We believe that (3 C and f3 c coincide, where f3 c is the critical inverse temperature for the dilute 
Ising model defined at (1.8). We also consider 

N = lim E${' f ' q + lim E$f w ' q 1 , (1.11) 

the set of inverse temperatures at which the infinite volume random media random cluster 
measure is not unique. It was shown in [Wou08], Theorem 2.3, that N is at most countable. 
Under the assumptions /3 > /3 C and (3 M, one can use a renormalization procedure (Theorem 
5.1 in [Wou08j) which gives a precise meaning to the notion of plus and minus phases and is 
hence a fundamental tool for the study of equilibrium phase coexistence. 



1.5 Surface tension 

Surface tension is another essential tool for the study of equilibrium phase coexistence. In the 
context of the dilute Ising model, it is a random quantity since it depends on the couplings J. 
We recall here some important definitions and results from [Wou09| . Let S"^ -1 be the set of unit 
vectors of R rf . Given n £ S^ 1 we let 

(d-l 

§> n = < ^^[—1/2, l/2]ttfc; (tii, • • • , u d-ii n ) is an orthonormal basis of R d 

U=i 
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(where ^ stands for the Minkowski addition) be the set of d — 1 dimensional hypercubes of 
side-length 1, centered at 0, orthogonal to n. Finally, given S G S n , x G R d and L,H > we 
denote 

H x , L , H (S,n) =x + LS + [-H,H]n (1.12) 

the rectangular parallelepiped centered at x, with basis x + L5 and extension 2/7 in the direction 
n. The discrete version of U is H = n n Z d and the inner discrete boundary of 1Z is 
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y G 72 : 3z G Z d \7e,z ~ y} . 



For any 1Z as in (1.12) we decompose &1Z into its upper and lower parts S + 72 = {y G 972. : 
(y — x) • n ^ 0} and d~TZ = {y G 972. : (y — x) ■ n < 0}. Then we call 

T>ti = |w G £1 : d + K & d~7e} (1.13) 
the event of disconnection between the upper and lower parts of dlZ, and 

•4 = -^be*& w (P*)- (1-14) 

the surface tension in 7Z. Surface tension is sub-additive and has a typical quenched value 

T ^ n) = J ^ r Wo,^N(5,n) in ^-Probability (1.15) 

that does not depend on 5 > nor on S 6 $„ (Theorem 1.3 in |Wou09j ) . It is positive for 
any f5 > f3 c (Proposition 1.5 in the same reference). We denote by J mm and J max the extremal 
values of the support of J. Since increases with J, surface tension can be as low as T™ n that 
corresponds to the constant couplings J = J mm , and as large as T^ ax when J = J max . According 
to the convergence in ( 1.15 ), r™ n and T^ ax also converge when 1Z = 1Zo,n,sn(<S , n) with N — > oo 
and we call their respective limits r mm (n) and r max (n). Surface tension can deviate from r|(n). 
Upper large deviations happen at a volume order (Theorem 1.4 in [Wou09]) and are irrelevant to 
surface phenomenon like phase coexistence. Lower deviations under r > r mm (n) occur according 
to the rate function 

W = ^-^^(r^ N>mi s, n) Kr), (1.16) 
see Theorem 1.6 in [Wou09]. The set 

Mi = j/3 > : 3n G S* -1 and r > such that Jjg jn (r|(n) - r) = o} (1.17) 
is at most countable, see Corollary 1.9 therein. 



1.6 Magnetization profiles 

In the following, C d stands for the Lebesgue measure on R rf and for the d—1 dimensional 

Hausdorff measure. The Z^-distance between two Borel measurable functions u, v : [0, l] d -> R 
is 

ll n — u IIl 1= / \u — v\dC d , 
J[o,i] d 
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and the set L 1 is 

ju : [0, l] d — > R Borel measurable, \\u\\ L i < oo| . 

In order that L 1 be a Banach space for the L 1 -norm, we identify u : [0, l] d — > IR with the class 
of functions {v : ||u — v\\ L i = 0} that coincide with u on a set of full measure. We also denote 
by V(u, 5) the neighborhood of radius 5 > in L 1 around kEL 1 . Given a Borel set U C 
we call 

'1 iixgU 
-1 ifsetf 



xv- x e~R d h+ 



the phase profile corresponding to U, and call V(U) the perimeter of [7 (as defined in Chap. 3 
of [AFPOO] ). The set of bounded variation profiles is 

BV = ju = xv ■ U C (0, l) d is a Borel set and V{U) < ooj . 

Bounded variations profiles u = xu G have a reduced boundary d*u and an outer normal 
n u : d*u —> S' 1 ' 1 with, in particular, Ti d ~ 1 (d*u) = V(U). As the outer normal n u defined on 
d*u is Borel measurable, we can consider integrals of the kind 

F q {u)= ! T q {n u x )dU d - l {x), VnGBV (1.18) 

J d*u 

that define the quenched surface energy of a given profile. When u = xu £ BV we also denote 
by F q {U) the surface energy of u. 

1.7 Initial configuration and gap in surface energy 



As described with further detail in the heuristics (Section 3.1), our strategy for controlling the 
Glauber dynamics is to start from some metastable initial configuration, from which a positive 
gap in surface energy must be overcome before the system can reach the pure plus phase. This 
metastable configuration is characterized, on the one side, by an initial phase profile no G BV, 
and on the second side, by a reduced surface tension r T on the boundary of uq. So a so-called 
initial configuration has actually two microscopic counterparts. First, to uq will correspond a 
set of initial spin configurations for the Glauber dynamics, while to the reduced surface tension 
r r we will associate a dilution event on the couplings J. 



Before we can define the set of initial configurations IC at (1.19), we need still a few more 
definitions. 

Definition 1.1. We say that a profile u = xu 1S regular if 

i. U is open and at positive distance from the boundary d[0, l] d of the unit cube, 

ii. dU is d — 1 rectifiable and 

iii. for small enough r > 0, [0, l] d \ (dU + -6(0, r)) has exactly two connected components. 

We recall that E C R d is a d — 1 rectifiable set if there exists a Lipschitz function mapping 
some bounded subset of onto E (Definition 3.2.14 in |Fed69j ) . It is the case in particular 
of the boundary of Wulff crystals (Theorem 3.2.35 in [Fed69| ) and of bounded polyhedral sets. 
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It follows from Proposition 3.62 in [AFPOO] that any u = xu regular belongs to BV and that 
dU = d*u up to a -negligible set. Finally, we call 

TP _ / ( u 0, T r ) G BV xC([0, l] d , R) : no is regular and there is 

~ \ e > such that, for all x G <9*n , r min (n£°) + e < r r (x) < t<?« ). 

Given an initial configuration (no,T r ) G IC, we define the reduced surface energy as 

F r {u)=f T r {x)dU d ' l {x)+ f T q (n x )ctH d - 1 (x) (1.20) 

which is obviously smaller than J- q (n). The reduced surface energy of the initial phase profile 
uq is 

.F>o)= / T r (x)dH d ~ 1 (x) (1.21) 

while the cost of dilution is 

F(u )= [ I nx {r r {x))dU d ~ l (x). (1.22) 

For any e > 0, we call C e («o) the set of sequences of phase profiles that evolve from uo to 1 
(pure plus phase) by jumps with L 1 norm less than e, that is 

r( \ \ t \ k G N ; v = u and v k = 1 \ , . 

C e (n ) = |^) i= o... fc • V < € {0, . . . , * - € BV and - ik|| x i < e / " ( ' 3) 

Finally we define the gap in surface energy or energy barrier for removing the droplet no, under 
the reduced surface tension r r : this is 

/C r (n ) = lim inf max F r {vi) - F r {u ). (1.24) 

£-5-0+ v£C e (u ) i 

1.8 The Glauber dynamics 

The Glauber dynamics is characterized by a family of transition rates c J (x,a) at which the 
configuration a changes to a x defined by 

y \ -o- y i£y = x. 

In other words, c J (x, a) is also the rate at which the spin at x flips. We make standard assump- 
tions on the transition rates, namely: 

Finite range There exists r < oo, the range of interaction, such that c J (x,a) is independent 
of a(y) when d(x, y) > r, and of J e when d(e, x) > r. 

Rates are bounded The rates are uniformly bounded from below and above: there are c m , cm G 
(0, oo) such that 

c m < c J (x, a) < c M , Vx £ Z d ,J £ J,a £ E. 




8 



Detailed balance The rates satisfy the detailed balance condition: for all a G £ and x G 
for all J G J, the product 



c J (x, a) x exp ( - ^ J{ x ,y}0- x a y j 



does not depend on a x . 
Translation invariance If, for some z E 7L d one has 

Jg = J z+e ,Ve E E(7L d ) and 4 = ^^x6^, 
then c J '(x, cr') = c J (x + z,a). 
Attractivity Given any a, a' G £ with cr ^ <t', the equality cj x = o' x implies 

a' x c J (x, a') ^ a x c J (x,a). 

Two important examples of the Glauber dynamics are the Metropolis dynamics, often used 
in computer simulations, for which 



\ y~x 



c (x,a) = max ^1, exp [ 2^ J{ x ,y)O x u y 

and the heat-bath dynamics 

c J (x,r) = n J {x} {a x = -r x \a y = T y , Vy ^ x) . 

Given the transition rates, one can proceed to a graphical construction of the dynamics, as 
follows: equip each site x G 7L d with a Poisson process valued on R + , with intensity cm (we 
recall that cm is a uniform bound on the rates of the dynamics). Consider now the time t 
growing from 0. When the Poisson process at x has a point at t, flip the spin at position x 
with probability c J (x,at)/cM- Because the flip rates are bounded, the determination of o~t(x) 
involves only finitely many sites and therefore the dynamics is well defined, even in the infinite 
domain 7L d . 

We call P J the law of the Markov process (a(t)) t ^o associated to this dynamics, and -P^(o) 
the law conditioned on the initial configuration <r(0). It is convenient to introduce the semi- group 
T J defined by 

[T J (t)f](a) = Ei(f(a(t)))- (1-25) 
The detailed balance condition makes the generator 

{L J f){a)= £ c J (x,a)(f(o-*)-f(o-)) (1.26) 

x&L d 

self-adjoint in L 2 (fi J ' + ), and ensures that the Gibbs measure fi J ' + is reversible for the dynam- 
ics. One way of quantifying the approach to equilibrium of the dynamics in infinite volume is 
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therefore the averaged autocorrelation 

A\t) = E([Va V , + (T J (t)vr )] A ) 



(1.27) 



E 



E 



[(T J (t)7T )(p)-^ + (a )] 2 d^ + (p) 



\T J (t)n -^(a )\\f 2{ ^ 



+) 



where 7ro : £ — > H is the function that, to the spin configuration a, associates ttq{(j) = ao, and 
A > is an arbitrary positive number. 

Although our work aims primary at describing the asymptotics of the averaged autocorre- 
lation, we also derive some upper bounds on the relaxation and mixing times in finite volume. 
To this aim, we introduce the Dirichlet form 

£ J A p (fJ) = W»a p {°Y j ^°)U{° x )- f (.<?)? 



(T,X 



and the spectral gap 



gapL 



J,A,p 



ml ( £ A^U1- f W ith VarJ'(/) + 



Varf(/) 

The inverse of the spectral gap is the relaxation time 

X J f' P = l/gapL^'". 
One fundamental property of the spectral gap is that, for any /, 

Vav J /(T^(t)f) < exp (-2t/7^ A '') Varf (/) 



(1.28) 



(1.29) 



(1.30) 



where T^' p is the semi-group corresponding to the Glauber dynamics restricted to A with bound- 
ary condition p. 

Finally, we define the total variation distance between two probability measures p and u, as 



\^ ~ HItv = sup /u(A) — zv(^4). 
a 



The mixing time is 



T 



J,A,p 



inf < t > : sup 

I ao:ao=p on A c 



j,p 
Ma 



TV 



^ e 



Given any function /, we have 



sup 

(Jo:<Jo=P on A c 



E^(f(a t ))-p J /(f)\ < exp(-[t/7^J) 



(1.31) 



(1.32) 
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2 Main Results 



2.1 Slow dynamics in infinite volume 



The main result of the paper is a rigorous lower bound on the averaged autocorrelation which 
confirm some of the claims of the foreseeing paper [HF87]. We recall that the set of (metastable) 



initial configurations IC is defined at (1.19), while the surface energy J 77 *, the cost of initial 



dilution Z r (uo) and the gap in surface energy JC r (uo) associated to (uo,t t ) S IC are defined at 



(1.20), (1.22) and (1.24), respectively. Now we define the exponent 

{u ,T r )&C:IC r (u )>0 /C r (u ) 



(2.1) 



for any A ^ 0. When no initial configuration (no, r T ) S IC leads to a positive surface energy gap 
1C(uq) > we adopt the convention that X\ = +oo. Under mild conditions the exponent X\ is 
finite and even bounded from above: 

Proposition 2.1. Assume that r min (n) < r*(n) for all TIGS' 1 - 1 . Then: 

1. For all (uo,T r ) G IC such that the boundary of uq is C 1 , the gap in surface energy JC r (uo) 
is strictly positive. 

2. If < P( J e = 0) < 1 — p c (d), for every A ^ there is C > such that, for any f3 large 
enough, 

C 
J' 



(2.2) 



Note that the above proposition is a corollary of Theorem 3.21 which is to be found, together 



with its proof, in Section 3.4.1 



Now we present our main result, which relates the decay of the autocorrelation of the infinite 
volume Glauber dynamics in the plus phase with the exponent X\ (we recall that M and A/j are 
defined respectively at ( |1.11[ ) and (1.17)). 

Theorem 2.2. For any j3 > j3 c such that ft ^ M Li Mi, for any A ^ 1, for any 5 > 0, for any 
t > large enough, 



A x (t) ^ t 



—X\—5 



(2.3) 



Remark 2.3. In general we have not been able to compute K7{uq) nor X\. Note however that 



we give another formulation of K7{uq) in Theorem 3.25 This alternate formulation is the key 
for computing /C r (uo) in two particular cases, see Sections 3.4.3 and 3.4.4 These computations 



also point out the fact that the constraint of continuous evolution in the definition of KJ{uq) at 



(1.24) makes, in some cases, the gap in surface energy bigger than if we only take into account 



the constraint of continuous evolution of the overall magnetization (see Lemma 3.29). 



Theorem 2.2 above was established some time ago during the PhD Thesis [Wou07j directed 
by Thierry Bodineau. Later on the author received the indications by Fabio Martinelli of how 
a corresponding upper bound could be established. The result is complementary to the former 
Theorem and confirms that the autocorrelation indeed decays as a power of time when d = 2. 
Furthermore, when < P(J e = 0) < 1 —p c (d), the exponent in both the upper and the lower 
bounds match up to a multiplicative constant as f3 — > oo. 
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Theorem 2.4. Assume that d = 2. For any A > 0, there is c > such that, for all (5 > j3 c , for 
all t > large enough, 



A x (t) < t 



-c/0 



(2.4) 



2.2 Slow dynamics in finite volume 

Our strategy for providing a lower bound on the autocorrelation also yields lower bounds on 
typical relaxation and mixing times in finite volume. Given d ^ 2, the law of interactions P and 
the inverse temperature /3, we call 

, IC r (u ) 
k = a sup — - — - 

(M ,T r )6lC 1 W 

= d/X 



We recall that the relaxation time T ve \ was defined at (1.29), while the mixing time T m i^ was 



defined at ( 1.31 ). 



Theorem 2.5. Assume d^2 and f3 > /3 C with @ ^J\fUj\fj. For any 6 > 0, 



MP r 



N 



-J,A? 
'rel 



> N' 



t-S 



limP (r m J ^' + > N«~ S ) = 1. 



It is instructive to compare this bound with the asymptotics of the relaxation time in the 
pure Ising model. When d = 2, T^ 2 ^' N ' + = N apart from logarithmic corrections, see [BM02J. 



For dilute models with < P(J e = 0) < 1 — p c (d), Proposition 2.1 states that k ^ c/3 for f3 large 



enough: this is another illustration of the fact that dilution makes the relaxation time much 
larger. 



3 Proofs 



3.1 Heuristics and organization of the proofs 



The object of Section |3.2| is the proof of the upper bound on the autocorrelation when the 
dimension is d = 2 (Theorem 2.4). The strategy for the proof, which we owe to Fabio Martinelli, 
relies on a uniform lower bound on the spectral gap in a square box with uniform plus boundary 



condition (cf. Theorem 6.4 in |Mar99| or (3.2) below). Some extra but classical work is then 
required to use that estimate for the infinite volume Glauber dynamics started from the plus 
phase. 



Section 3.3 concentrates the most complex part of the paper and is dedicated to the proofs 



of Theorems 2.2 and 2.5 The corresponding heuristics are derived from [HF87] . where it was 
already suggested that dilution and reduction of surface tension, which has a surface cost, could 
trigger metastability for initial spin configurations corresponding to the minus phase in the 
region surrounded by the diluted surface (i.e., with our notation, the region where Uq = — 1, up 
to an appropriate scaling). In the present work, we have formalized that idea of reducing the 
surface tension along an initial contour with the set of initial configurations IC defined at ( |1.19 ) . 
The concept of energy barrier has lead to the definition of the gap in surface energy K7{uq) at 
(1.24). The reader will note that we have introduced, in that concept of energy barrier, the 
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requirement that evolution of the phase profile is (almost) continuous in the L 1 norm (with 
initial configuration no and final configuration 1, corresponding to uniform plus phase). Once 
these concepts are defined, remains a substantial work for relating the concepts and definitions 
to the physical phenomena. This is done as follows. First, in Sections |3 .3 . 1 1 and 3.3.2, we define a 
microscopic counterpart for the (uo,t t ) dilution, the event G((TZi)i=i.. n ), which depends only on 
the interaction strength J e . We prove that this event has the expected probability (the cost for 
dilution) and also we establish, conditionally on G((TZi)i=i.. n ), upper and lower bounds on phase 
coexistence. In Section 3.3.3 we introduce a certain bottleneck set, that corresponds with the 
phase profiles with the highest cost along all relaxation paths from uo to 1. Then we express the 
probability, under the equilibrium measure, that the actual phase profile is in the bottleneck set 
in terms of its reduced surface energy, and in turn we relate the infimum of the reduced surface 
energy in the bottleneck set to the energy barrier JC r (uo). In Section 3.3.4 we relate the former 



estimates on the equilibrium measure to the Glauber dynamics and establish an important 



result on the dynamics in a finite box (namely, Proposition 3.18). That Proposition states that, 



conditionally on the event of dilution, the average magnetization remains significantly different 
from the plus phase magnetization (i.e. a droplet of minus phase remains), with a probability 
corresponding to the reduced surface energy, until a time determined by the energy barrier. We 
conclude the proof of Theorem 2.2 and Theorem 2.5 in Section 3.3.5 For the proof of the first 



Theorem we relate the autocorrelation to the evolution of the magnetization in a finite volume. 
For the proof of the second Theorem, we decompose the box of interest into many smaller boxes, 
so that the dilution event must occur in one of these boxes with a probability close to one. 

Finally, Section [3. 4| is concerned with geometrical estimates. By decomposing the difference 
between the initial droplet and the current droplet into many small droplets, we are able to 
use the assumption that the initial droplet has a smooth boundary and prove that the energy 
barrier is strictly positive. We also prove that, under appropriate assumptions, the exponent X\ 



is bounded by C/f3 (Proposition 2.1). Then in Section 3.4.2 we derive another formulation of 



the energy barrier, which shows that the energy barrier can be computed assuming not only the 
L 1 continuity of the droplet evolution with time, but also the 'K d ~ x continuity of the part of the 
droplet surface that corresponds with the initial droplet contour (where dilution takes place). We 
establish that alternative and informative formulation by providing an interpolation between any 
two phase profiles that is continuous for both volume and surface measures, with the additional 
property that along the interpolation, the surface energy does not exceed the maximum of the 



surface energy among the two interpolated phase profiles. Finally, in Sections 3.4.3 and 3.4.4 



we use that alternative formulation to compute the energy barrier in specific cases. These 
computations put in evidence the fact that the energy barrier, as defined by /C r (ito), is higher 
than the energy barrier computed with a constraint on the overall magnetization. 



3.2 Upper bound on the autocorrelation 



Here we give the proof of Theorem 2.4 As said above, the scheme of proof was suggested by 
Fabio Martinelli. We call 



m 



(7 



Tl' + (t)n )(a) 



J ' + (^o) 
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therefore // (a) is the mean value of the spin at the origin under the Glauber dynamics performed 
until time t, initiated with the configuration a, while is the corresponding quantity for 

the dynamics restricted to A, with plus boundary condition on A c . Using these notations we 
can write the averaged autocorrelation A x (t) defined at (1.27) as 



A x (t) = E [Va V , + (//)] 



Lemma 3.1. One has 



(3.1) 



Va V , + (//) < //(+) 



m 



where + means the constant plus spin configuration. 

Proof. The attractivity of the dynamics implies that //(<r) ^ //(+) f° r an Y Then: 



Va V , + (//) 



f{(*)>ft J (.v) 



{ft (*) ~ ft (V)) d/-+ (^d^+fa) 



(// (+) " // fa)) <W'+ (r?) 



(x 2 s£2x,Vxe[0,2]) 
attractivity 



And the last term equals //(+) — w^. 

We apply then standard controls on the spectral gap: 



(//M- //(77))d//' + (a)d/' + (t?) 



□ 



Lemma 3.2. There exist positive and finite constants G\,C2,C$ such that, for all J G ', aZZ 
i > and /3 > 0, 



/t, t (+)-/L,H < e 



C 1 t3N d -C 3 tN- d e- c 2fiN d 1 



Proof. The difference /f (+) — // ( 



according to (3.15) in |Mar99j does not exceed 



is not larger than twice the ||.||oo norm of // , which 

A N ,t 



J A N ,t 



inf /j J ^ + (a) exp ( —t gap(L 



-1/2 



Fo| 



The conclusion comes then from the general lower bound 

gap(L J ' A ^+) ^ c m (2iV + l)- rf e - C7 ^ Jvti - 1 ; 
cf. Theorem 6.4 in |Mar99j . 

An easy consequence of monotonicity is the next Lemma: 



(3.2) 

□ 
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Lemma 3.3. For any J € J ', one has 

^ 2/'+(CX) 

where C\ is the event that there exists a contour of plus spins in A, around the origin. 



Proof. The inequality f\ t {—) ^ /x^' + (o"o) is immediate. Then we remark that, conditionally on 
C\, the expectation of o~o under fi J is at least ^ + (00 )• Otherwise it is at least —1. Hence 

m J > /' + (C A )x^ + (ao) + / u J ' + (CX)x(-l) 
^ ^ (ao) - 2/'+ {C%) . 

□ 

The average of n J ' + [C\ ) decays exponentially fast with the size of A under the slab perco- 
lation assumption when d = 2 (see (1.10) for the definition of j3 c ). When f3 is large, we can also 
give quantitative estimates. 

Lemma 3.4. Assume d = 2 and (3 > (3 C , then there is c > such that, for N large enough, 

E//>+(qJ < exp(-ciV). (3.3) 

Proof. We use the FK representation of the spin model and the renormalization framework of 
|Wou08| . The event is realized with conditional probability one if both of the following 
occur in the edge configuration: 

i. There is an infinite w-open path issued from A N / 2 

ii. There is a surface that lies in Ajv \ A N / 2 f° r which all points are w-connected. 

Since d = 2, the second point reduces to finding a circuit of open edges inside An \ A N / 2 . We 
can cover A^v \ An/ 2 with 16 blocks of side-length N/2. When all these blocks are good in the 
sense of Theorem 2.1 in [Wou08| . the second point is realized, and this occurs with a probability 
1 — exp(— cN) as we assume that (3 > f3 c . Similarly, the first point is realized as well with a 
probability 1 — exp(— cN) and (3.3) follows. □ 

Remark 3.5. In larger dimensions, under the assumption that P(J e ^ e) = 1 for some e > 0, 
one can use Peierls estimates to prove that ^C? ^ ^ exp(— cj3N) for some c > 0, for 

any f3 large. Still, this is not useful for generalizing Theorem |2.4| as can be seen in the final 
optimization below. 



Now we give the proof of Theorem 2.4 



3.1 



Proof. According to Lemmas 
monotonicity of the dynamics, we 



3.2 
have' 



3.3 



and to the inequality //(+) ^ f^ t (+) due to the 



Va V , + (//) ^ 2e C lf 3N d -C 3 tN-*e-W ~ +4 /.+ (q 
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We recall the assumption that d = 2 and j3 > f3 c . According to Lemma 3.4 there is c > such 
that, for all ./V large enough: 

E//'+(qJ ^ exp(-ciV), 



hence Markov's inequality implies that 



P \ji J >+ (C|J ^ exp(-o/V/2)J < exp(-edV/2). 
Now we define iV by the relation 

C 2 pN d ~ l = (l-tf)logt. 

For large enough t, we have 

Vav >+ (//) < 4(exp(-c7V/2) + /, J ' + (q 



Combining this with (3.1) and ( |3.4[ ) gives 

A x (t) < 8 A [exp(-cA./V/2) + exp(-dV/2)] 



(3.4) 
(3.5) 



which, according to the assumption that d = 2 and to the definition (3.5) of N, establishes the 
claim dOb. □ 



3.3 Lower bound on the autocorrelation 



The subject of the present Section is the proof of Theorems 2.2 and 2.5 It is organized as 
follows. In Section 3.3.1, we define the notion of covering of the border of a magnetization 
profile by rectangles. Then, in Section 3.3.2 we define the event of dilution according to some 
initial profile (uq, t t ) G IC and prove that it has the expected probability. We show how previous 
result fro m [Wou 09] apply for the probability of phase coexistence, given the event of dilution. 
In Section 



3.3.3 



we show that phase profiles evolve continuously in L 1 and relate this property 
to the bottleneck and to the gap in free energy /C r (uo). Finally in 3.3.5 we conclude the proof 
of Theorems 12.21 and 12.51 



3.3.1 Covering of the boundary of phase profiles 

As in the work [Wou09j, the coverings of the boundary of macroscopic phase profile are a 
fundamental tool for relating the macroscopic shape of the magnetization to the microscopic 
spin system. The definition that we present here is more restrictive than in |Wou09] since it 
takes into account the set IC, together with a new parameter 7 (fourth line in (v) below). 

Definition 3.6. Let (tto = XUoi tV ) £ IC and u G BV, together with 5, 7 > 0. We say that a 
rectangular parallelepiped 1Z C [0, l] d is (uq, r r , 5, 7)-adapted to u at x 6 d*u if: 

i. If n = n" is the outer normal to u at x, there are S £ S n and h G (0,6] such that, either 
TZ C (0, l) d (we say that 1Z is interior) and 

K = x + hS + [±Sh]n, 

either TZn d[0, l] d 7^ (we say that 1Z is on the border) and x G <9[0, l] d , n is also the outer 
normal to [0, l] d at x and 

n = x + hS + [-5h,0]n. 
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ii. We have 



U d - X (d*undTZ) = 0, 



1 



1 - t^t^ 1 (d*unn) 



h d 



and 



iii. If x : IR d 
namely 



then 



^ 5. 



h Jd*umz 

{±1} is the characteristic function of the half-space above the center of ft, 
+1 if (z - x) ■ n > 



X{z) 



-1 else 



,Vz E R a , 



1 



25h d 



\ X -u\dH d < 5. 



7v 



iv. If ft n <9[0, l] d / 0, then ft does not intersect dU . 

v. If ft intersects dUo, then x E dUo and 

1 -/d-i 



jTT" 1 ((d*uAd*« ) n ft) 



r x 



1 



T r dH 



d*u nll 



and 



J»(t p (*)) 
4(r r W-7) 



1 



1 



d*u mi 



^o(r^(.)- 7 )^- 1 (z) 



d*u nn 



«c 5. 



vi. If ft intersects 9C/o, then the enlarged volume 

ft' = ft + 5(0, 2\/d/i 2 ) = {z£E J : d(z, ft) 2\/d/i 2 } 

satisfies 



T w < *- 1 (d*u n ft' \ ft) < <5. 



(3.6) 



Note that conditions i to iii above mean exactly that ft is <5-adapted to d*u at x E d*« in 
the sense of Definition 3.1 of [Wou09j. Also, by A we mean the symmetric difference. 

Definition 3.7. Let (ito, r r ) E IC and u E BV, 5, 7 > 0. A finite sequence (ftj)i=i... n of disjoint 
rectangular parallelepipeds included in [0, l] d is said to be a (uq, r r , 5, 7)-covering for d*u if each 
fti is (no, T r , 5, 7)-adapted to u and if 



u A ~ x a*u \ Q fti ^ 5. 



(3.7) 



i=i 
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The main result of this Section is: 

Proposition 3.8. Let {uq,t t ) G IC and u G BV, together with j,5 > 0. There exists a 
(no, T r , <5, ^-covering for d*u. 

Proof. As the proof follows a classical argument, we will only give here the main steps in the 
proof. First, we claim that, for any u, uq G BV, for Ti^ 1 almost all x G d*u fl d*uo, 

lim -4t n d ~ 1 ((d*uAd*u )nB(x,r)) = 0. 

r-s>0+ r a 1 

This can be proven by applying, for instance, the Besicovitch derivation Theorem (Theorem 
2.22 in |AFP00| ) to the Borel measurable function 



f :x£~R d ^ 



1 if x G d*uo 

2 else, 



with the result that for % d -almost all x G d*u n d*uo, 



lim l —r— [ fdU d - x = fix) 

r^0+ Ol d _ir a Jd* ur] B(x,r) 

where a d _i = H d ~ x ({x G 5(0, 1) : x ■ e d = 0}) . Therefore, 



lim 

r->0+ 



(3*u n B(je, r)) ?^ \ 9*n ) n 5(x, r)) 



As the first term goes to 1 already as r — > 0, we conclude that 

1 



lim - — t Zd z I n d - 1 ((d*u \ d*u ) n S(x, r)) = 



for H d 1 - almost all x G d*u Pi <9*uo- The claim follows as it and ito play a symmetric role. 

Second, we remark that more generally, all the absolute values in point (v) of Definition 



3.7, and also left-hand side of inequality (3.6), have a zero limit as h — > for % d 1 -almost all 
x G d*u n dUo. This is a consequence of the strong form of the Besicovitch derivation theorem 
(Theorem 5.52 in [AFPOOj ^. 

The two above facts enable, like in the proof of Theorem 3.4 in [Wou09|, the use of the Vitali 
covering Theorem (cf. Definition 13.2 and Theorem 13.3 in |Cer06j ) to conclude to the existence 
of the desired covering. □ 

3.3.2 The event of dilution 

The purpose of this Section is to define the event of (microscopic) dilution and to prove several 
important properties about phase coexistence, given the event of dilution. We recall a notation 
from |Wou09| : when 1Z is a rectangle and N > 0, we let 

U N = NTl + z N (K) 

where zn(1Z) G (—1/2, l/2] d is such that the center of 1Z N belongs to Z d . 
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Definition 3.9. Let (uq, t t ) E IC, 5, 7 > and let (TZi)i=i... n be a (uq, r r , 5, 7)-covering of 8*uq. 
The event of dilution on this covering is 



i )i=l...n ) 



J E J :t£„ < r r (xi),Vi = 1 . . . n} 



(3. 



The reader will check that the event of dilution affects only the random variables J e with e 
at distance at most N5 from the boundary Nd*uo- 

The first property of the event of dilution is that it happens at the expected rate. Namely, 

Lemma 3.10. Let (uo,T r ) E IC and £ > 0. For 5 > small enough, for 7 > arbitrary, if 
(7£j)j=i... n is a (uo,T r , 5, ~f)- covering ofd*uo, then 



P (0 ((fti)i=l...n)) > exp 
/or any iV iarge enough. 



-N 



d-l 



Proof. We recall that X r («o) was defined at (1.22), see also (1.16). Since the TZf are disjoint 
for large enough N, we have 



liminf-^logP(S((7^ =1 . 



n ^ 

■»)) = ^Ejv^T^fe^^)) 

n 

= -E^«*(^)) 



in view of (1.16). Note the role played by the assumption r r (x) > r min (n"°) in the definition 
of IC at (1.19). The properties of the covering (Definition 3.7 and point (iii) in Definition 3.6) 
imply the claim for 5 > small enough. □ 

Then we show that the dilution has the expected impact on the probability for phase coex- 
istence. This is expressed with the two complementary Propositions |3 . lT] and 3.12 We have to 
recall here the definition of the magnetization profile 



M 



K 



[o,i] c 

x 



K 



73 S 



■1,1] 



zeAjvDAi 



Or) 



0~z 



(3.9) 



where K E N* is the mesoscopic scale and 
i(x) 





'Nxi 




'Nx d ~ 




( 


K 


j . . . , 


K 


) 



and Ai = Ki + {1, . . . , K} c 



(3.10) 



Hence, unless x is too close to the border of [0, l] d , Mk{x) is the magnetization in a block of 
side-length K that contains Nx. Theorem 5.7 in [Wou08] provides a strong stochastic control 
on Mk when j3 > j3 c . In particular, when K is large enough, at every x the probability that 
M.k(%) is close to either m$ or — rap is close to one under the averaged measure E/z^. The 
event Ai^c/mp E V(uo,e) means therefore that the system is close to plus (resp. minus) phase 
at Nx when uq(x) = 1 (resp. uq(x) = —1). 
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Proposition 3.11. Assume that (3 > (3 C and (3 ^ M. Let (no, r r ) G IC and e, £ > 0. Ze£ 5, 7 > 

and consider (J^i)i=\... n o> (uo,T r ,5, 7)- covering of 8*uq. Then, if S > is small enough (and 7 
arbitrary), for K large enough, 



lim P I 



M 



K 



G V(«o, e) ) ^ exp ( —N d ~ l (T r (u ) + £) 



G {(Hi) 



i=l. ..n 



1. 



Proof. We use the notations of Proposition 3.9 in |Wou09| . We recall that D^' 5 is the event of 



w-disconnection around Nd*uo and that £?t' S is the set of edges close to Nd*UQ. We let then 



J,w,+ 



F J N = inf * A 



M 



K 



mp 



G V(u ,e) 



uj = -/r on £, 



N,S 
U 



Proposition 3.9 in |Wou09| and the definition of the covering imply 
i. For 5 > small enough, 



lim inf 



inf 



jv^oo Jeg((Hi)i=i...„) N d 
ii. For 5 > small enough, for K large enough 



lim P|Fjv< 



1 



0. 



(3.11) 



(3.12) 



The definition of FjL yields 



J,+ 



M 



K 



mp 



6V( U0 , £ ))^^;(^) 



hence, using (3.11) we obtain: for large enough N, 



i)i=X...n, 

> V(F J N > l\g{(TZi) l=1 ... n )) 



Yet, the variable is independent of the J e with e G £jjf ■ Thus it is as well independent of 
the dilution Q ((11, 



N 

i )i=l...n ) 



-U 

and (|3.12l) yields the conclusion. 



We also recall from [Wou09j the definition of the compact set 

BV a = {u = xu G BV : V(U) ^ a} 
where V(U) is the perimeter of U. 



□ 



(3.13) 



Proposition 3.12. Assume f3 > j3 c and (3 AAuW}. Let (u ,T r ) G IC, a > 0, £ > 0. There 
exists 7 > such that, for any u G BV a there is e > such that, for any 6q > small enough, 
for any (uo,r r ,5q, ^-covering (TZi)i=i... n of 8*uq, for any K large enough, 



lim P(^ + (— eV( U , £ )| ( ' X P ( - A'"'~ J (F'( a.) — ( ) 

N^-oo \ N \ mp 



g((n 



i )i=l...n, 



1. 
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The proof of Proposition |3.12 is based on Proposition 3.11 in [Wou09j that relates the 
probability //^' + (A4k /mp S V(u,e)) to the Z^-notion of surface tension 



~J,5,K 
r N1Z 



N1Z 



NK 



M 



K 



X 



^ 28C d (K) 



(3.14) 



together with Propositions 3.12 and 3.13 in [Wou09| that compare the tw o defi nitions of surface 

-J,S,K 



tensions and Before we complete the proof of Proposition 



typical value of f^^ K given the event of dilution. 



3.12 



we examine the 



Lemma 3.13. Assume (3 > f3 c and (3 jVu A/j. Let (uo,T r ) G IC and 7 > 0. For 5 > small 
enough, the following holds: for any u G BV and any 1Z that is (uq,t t ,5^)- adapted to d*u at 
x G d*uC\d*uo, for any (uo,r r , <5o, 7 ) -covering (TZi)i=i... n ofd*UQ such that So G (0, h 2 ), for any 
K large enough, then 

^ — Lj logP (fjgf < T r (x) - 4j -y\Q m) i=1 ... n )) < 

where c' dS = c^S + c^s is the sum of the constants that appear in Propositions 3.12 and 3.13 in 
lWou0 9l. and is arbitrary small as S — > 0. 



Note that h in the statement of Lemma 3.13 refers to the largest dimension of 7Z, see (i) in 
Definition 13.61 



h 



Vdh 2 




Figure 1: The scale of dilution 5q. 



Proof. (Lemma 3.13). We consider (7£j)j=i... n a (uq, r r , <5q, 7)-covering for d*UQ. Thanks to the 



product structure of P, for large enough iV the surface tension f£f} K is independent of the r. 



such that TZi n 1Z = 0. Hence, for large enough N, the conditional probability 

P ( f NK K < ^0) ~ c 'd,8 - 7 Q ((ftt)i=l...fi)) 

equals 

P {jNK < r r (x) - c' d>s - 7 | < r p (zi), Vi : 7^ n K £ <b) , 
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which is not larger than 



P If 



PN 



t J,6,K 
Nil 



< T T {x) -C' d>s -7) 



P (r^v < r r (xi),Vi iKiMl^ 



According to the definition of Cj { , to Propositions 3.12 and 3.13 in [Wou09j (here we use the 

assumption that f3 > f3 c and j3 A/*) and to the definition of / n at (1.16), for large enough K 
we have 



A? 



(3.15) 



We will show that the right-hand side of (3.15) is negative, for small enough 5 > and 5o ^ h . 
Thanks to (vi) in Definition 3.6, one sees (Figure [TJ that the right-hand side of (3.15) is not 
larger, for S ^ 1/2 and 5q ^ h' 2 , than 



[/„, (r r (*)) - I nz (r r (z) - 7 )] dU d -\z) + M^' 



,<f-l 



where 



M = 3 + 2 sup I^o (r r (x)) 



is finite thanks to the definition of IC at (|1.19|). Now we give an upper bound on the integral. 



Since I n is convex, its slope is non-increasing and therefore 
[InAr r {z))-I nz {T r {z)- 1 )]dH d -\z) < 



d*u mz 



d*u mz 

d-l 



(l/2)h° 



inf / n (r 9 (n) - 7) 



Now we show that inf^ggd-i I n (r q (n) — 7) > 0. If not, we can extract a converging sequence 
rik — > n in the compact set with I nk (r 9 (n^) — 7) — > as A; — > 00. The Fenchel-Legendre 
transform r A (n) = inf T>T Diw ra \{AT + 7 n (r)} satisfies therefore 



r A (n fc ) < Ar 9 (n fc ) - A7 + o (1), VA > 0. 



k— >oo 



Since r A and r q are continuous (Proposition 2.5 in [Wou09j), we obtain in the limit r A (n) ^ 
Ar 9 (n) — A7, and therefore, by duality of the Fenchel-Legendre transform I n (r 9 (n) — 7) = 
0, which contradicts the assumption f3 ^ Mi- The claim follows for any 5 < 1/2 with 5 < 
mf neSd - 1 I n ^(n)- 1 )/(2M). □ 

Proof. (Proposition 3.12[ ). Let 7 < £/(2a). We take 5 > small enough so that we can use 
with c' dS ^ 7. Let u G BV a and consider (7lf)i=i...n(u) a ( n o> t t , S, 7)-covering of 



Lemma 



3.13 



d*u. Proposition 3.11 in [Wou09j states that 



N d 



n(u) 



u\d-l ~J,5,K 
T NKf 



i=l 
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for e > small enough. According to the definitions of 7 and 5 and to the properties of the 
coverings, provided 5 was chosen small enough (still independently of u E BV a ) we have 

n(u) 
i=l 



Finally when we take 5q < mm^^(hY) 2 and K large enough, Lemma 



3.13 



concludes the proof. 

□ 



Finally we remark that dilution has little influence on the overall magnetization 



(3.16) 



Proposition 3.14. Assume that (3 > /3 C and (3 ^ M. Let (uo,r r ) E IC and e > 0. Let 5, 7 > 

and consider (TZi)i=i... n a (uq,t t , 5, 7) -covering of 8*uq. Let N > and 



£ N = e(((]k? ) (J (A* + (Z d \ {0} 



N + 



the set of edges in some Tif , or in some of the translates of An by z(N+ [vJV]), for z E Z \{0}. 
Then, if 5 > is small enough, 



limP (/i J ' + ("t-Ajv) ^ "fy? — z\Je = 0, Ve E £at) 



1. 



Note that fj, J,+ (mA N ) increases with every J e , therefore the condition that J e = 0, Ve E £n 
is the worse condition that one can consider. 

Proof. The proof is based on the renormalization procedure established in Theorem 5.10 in 
|Wou08| . As it is similar to that of Proposition 3.9 in [Wou09j we only sketch the argument. We 
cover E(1j d )\£]\[ with blocks with side-length Ljy = [V^V]- To these blocks is associated a family 
of independent random variables (f{ E { — 1,0,1} (independent of the J e ,e E £n), the local e- 
phase, that is zero with probability less than e~ c ^ . Write /j, J,+ (mA N ) = limj^ /i J ~' + (rriA N ) and 
take some finite M. A simple Peierls estimate shows that, given the condition J e = 0, Ve E <?jv, 
the e-local phase associated to any block in A^ is +1 with a probability going to one as N — > 00, 
uniformly in N. As described in Theorem 5.10 in |Wou08| . this event implies that 



^ (mp-e) (l-^Vol(^) 



i=i 



Therefore we just need to take e > small enough (named 5 in Theorem 5.10 in [Wou08j) and 
5 > small enough so that Ya=1 Vol(7£j) ~ 5'H d ^ 1 (d*U( i ) is negligible. □ 
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3.3.3 The bottleneck 



Here we focus on the bottleneck in the dynamics. We recall that the set C e (uo), the set of 
sequences of profiles that evolve from uq to 1 with jumps in L 1 -norm less than e, was introduced 



at (1.23). Given v = («i)i=o...fc £ C e (uo), we call 



nun <; i : inf J~ r (u) = max inf J~ r (u) 



argmax £ (v) 

Now, given (uq,t v ) E IC, we define the e-bottleneck set as 



[J y( u argmax 2e »> e )- 
v&C E (u ) 



(3.17) 



When (no, r r ) are clear from the context, we simply write B e for B s (uq, t t '). Note that our motiva- 
tion for the above definition is that the e-enlargement of B e is V(B e , s) = \J v <=c e (u ) ^( t, argmax 2e («) > 2e), 



a fact that helps in the proof of Lemma 3.17 below 



Now we state three Lemmas related to the bottleneck set. Lemma 3.15 gives the asymptotics 
of the probability that the phase profile M.K/mp belongs to B £ . In Lemma 3.16 we show that 



T r is lower semi-continuous, a requisite for the proof of Lemma 3.17 in which we show how the 



gap in surface energy K7{uq) defined at (1.24) is related to the e-enlargement of B £ {uq,t 



Lemma 3.15. Assume f3 > /3 C and (3 ^ jVuA/}. Let (uq,t t ) G IC and £, s > 0. There exists 
7 > such that, for any (uQ,T r , 5, 7) -covering (Jti)i=i,„ n of d*u$ with 5 > small enough, for 
any K large enough, 



-N*- 1 ( inf F - t 



i i=l. 



.n) 



1. 



Proof. The exponential tightness property (Proposition 3.15 in |Wou09j ) tells that there exists 
C > such that, for every e' > 0, for any K large enough, 



hmsup log E/4>; V(BV a , e>] 



-Ca 



(the set BV a was defined at (3.13)). An immediate application of Markov's inequality shows 
that 



1 



N d 



' '^V(BV a ,< 

"2-/3 



Ca 
~2~ 



with a probability at least 1— exp(— CaN /3), for large iV. Since the surface cost oiQ{{1Zi) 
is bounded, we have therefore, for any e' > 0, for any a large enough, for K large enough: 

Ca 



i=l. ..n J 



]miI> { Ma'J (^^ V ( BV - e ')) ^exp 



-iV 



<?((7^) 



i=l. ..n. 



1. (3.18) 



Now we take some £, e > and a > large. We take for 7 > the one given by Proposition 



3.12 (it does not depend on e), and for any u G BV a we denote by e^(u) > the parameter e 
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given by the same Proposition, and e(u) = mm(e, e^(u)). Since the set BV a is compact, it is 
covered with a finite number of balls 



BV a C \JV( Ui ,£(ui)) 



i=i 



where U{ € BV a , for all i £ {1, . . . , re}. Since the right-hand set is open while BV a is compact, 
there is e' > such that 



V(BV a ,e') c\JV( Ui ,e( Ui )) 



i=i 



and we can write 



B e nV(BV a ,eO C |J V(u u e( Ui )) 

i^n:Ui<=V(B £ ,e) 



It follows then from Proposition 3.12 that, for any 5 > small enough, for any (uq, t v , 5, 7)- 
covering (TZi)i=i... n of 8*uq and any K large enough, 



limP 

N 



Ma 



N \ m/3 



eB e r\V(BV a ,e')) sC 



nexp (-N d 1 (minj <n:u . eV(Beje ) J 7- ^) - £)) 



0((fti)i=l...n) = 1. (3.19) 



Combining (3.19) with ( |3. 18 ) proves the Lemma, provided that a > was chosen large enough. 

□ 



Lemma 3.16. For any (uq,t t ) £ IC, the functional T r defined at (1.20) is lower semi- 
continuous. 

Proof. We show the lower semi-continuity as an application of the covering Proposition (Propo- 
sition 3.8). Let u E BV and 5, 7 > 0, and consider a (uq, r r , 5, 7)-covering (1Zi)i=i... n for d*u. 



Since the IZi are disjoint, for any v £ BV we have 



•^0) > J2 f. T r (x)dU d - l {x)+ ! 

i=i J'R- l nd*u nd*v J(ii % 



\d*uo)nd*v 



> y U A - X (Kind*v) inf r r + Y [ T q {n v x )dU d -\x 



i-.TlindUo^ 



i:^ 4 n9(7o=0 ' 



since r r (x) ^ T q (n x ). Thanks to the lower semi-continuity of the surface energy in open sets 
(Chapter 14 in |Cer06| ). the quantities 'H d ~ 1 ('R,ind*v) and f-^. n Q* v T q (ri" x )d'H d ~ l (x) become not 
smaller than their value at u when v converges to u in L 1 norm. Hence: 



lim inf r(.) > £ {l-5)h d ~ l inir r + £ (1 - 



n, 



which is arbitrary close to J- r (u) for small 5, thanks to the uniform continuity of r r . 
Lemma 3.17. For any (uq,t t ) £ IC ; 

K7{uq) = lim inf J 71 ' (u) — J 71 " (uq) . 

£->0+ ueV(B e (u ,T r ),e) 



□ 



(3.20) 
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Proof. We prove first that K7(uq) is larger or equal to the right-hand side. Let v E C £ (uo). 
Then, 

max J^Oui) ^ max inf T r 

i i V(v it 2e) 

^ inf F r (u) 

and, when we optimize over v E C £ (uq) and let e — >• we obtain the inequality 

/C r (u ) + -F r (n ) ^ lim inf J" r (n). 

e^0+ MeV(B s (M ,T r ),e) 

Now we prove the opposite inequality. Given any e > 0, there is n G V(i3 e (uo, r r ), e) such that 

r» < inf F(w) + s. (3.21) 

weV(B e (Mo,r r ) i e) 



According to the definition of B e (uQ,T r ) at (3.17), there is v = (^i)j=o...fc in C e (uq), that is an 
e-continuous evolution such that vq = uq, Vk = 1, that satisfies 



•7^ (it) #s max inf J^Cw). 

i=0 tueV(ui,2e) 

Now, for each i = 0, . . . , k we consider t^ +1 G V(t>j, 2e) such that 

r{v' l+l ) < inf ^H + e. 

tueV(uj,2E) 



(3.22) 



(3.23) 



We let also v' = no and ujj.+2 = !• Clearly, the evolution ?/ = (f-)i=o...fc+2 is 5e-continuous. 
Therefore, we have 

inf maxJ^f/) < max 7^ («,■)• 

«"eC 5e («o) < i=o 

The maximum in the right-hand side does not occur at i = k + 2 since J-" r (l) = 0. We also have 
the bound 

F{v> G )-F{v> x ) < J"" (no)- mf o ^ r H 

therefore, 



fc+1 



inf maxJ^^ ) ^ maxj rr (f i ) + ^(uq) — inf ^(w) 



v"GC 5E (u ) i 



i=l 



™eV(uo,2£) 



^ max inf J 77 " (w) + e + ( ^(uq) — inf F r (w 



i=0 u)eV(«i,2e) 



weV(«o,2e) 



< -F'(u) + e + J^ r (n ) - inf F r (w) 

X w&V(u ,2e) 



inf F(w) + 2e + J^(no) - inf .F7n/ 

iu6V(5 e (w ,T r ),e) V toeV(u ,2£) 



where the second line is due to the definition of v' i+l at (3.23), the third one to (3.22) and the 
last one to (3.21 ). The lower semi-continuity of J- r (Lemma 3.16) imply that the last term goes 
to as e — > 0. Therefore taking e — > ends the proof. □ 
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3.3.4 Intermediate formulation of the metastability 

The aim of this Section is the proof of an intermediate and useful formulation of the metastability: 

Proposition 3.18. Assume ft > j3 c and (3 M L)J\fj. Let (uo,r r ) E IC, £ > and e > small 
enough. Then, there exists 7 > such that, for any 5 > small enough, for any (uo , t v , 5, 7) - 
covering (Tli)i=x... n ofd*u , 



N 

d-l 



i )i=l...n, 



1 



where K N = exp(iV d - i ()C r (u ) - £)) 



Important keys for the proof of Proposition 3.18 are Lemmas 3.15 and 3.17 about the bottle- 
neck of the dynamics, together with Lemma 3.19 below on the so-to-say "continuous evolution" 
of the magnetization profile. 

Lemma 3.19. Let (3 > and e > 0. There is c > such that, for any N large enough, 



sup S npP J ^' + (\\MK(<7(t))-M K no))\\ Ll >e) < exp 

te[0 >£ /(4c M )] J 

where cm is a uniform upper bound on the rates of the Glauber dynamics. 



-cN c 



Proof. (Lemma |3.19|) . The L 1 distance \\M K (cr(t)) - M K (<r(0))\\ L i is bounded by 2/N d times 
the number of jumps of the Glauber dynamics. These jumps occur at a rate bounded by cm- 
Therefore, 



P J ^ + (\\M K (a(t))-M K (a(0))\\ Ll 



> e) sC Pi X> 



eN c 



where X is a Poisson variable with parameter CMtN d ^ eiV d /4. Cramer's Theorem imply the 



claim. 



□ 



Proof. (Proposition 3.18). In view of Proposition 3.11 it suffices to prove that 



lim inf Pip ^ — e 

N t^K N \ 2 



9 ((Wi)i=l...n) 



where 



j ;+ / Mx/nip G V(uq,s) and 



Pa\ 



r £"(*) m A JV > rnp - 2e 



So we focus on the proof of (3.24). We remark that, for any initial configuration p, 
P J p An ' + (m AN (a(t)) ^mp-3e)<e (1%+ (t) m Ajv )(p) < m p - 2e 

since (T^ (t) m Ajv )(p) = E J p ^ N ^ + (m Ajv (cr(t))) and m AN (a) < 1, for any a. Therefore 
J 



P 



j. + ( M.K/mp G V(uo,e) and 



^ 1 p^Ajv 

*5 ' J,+ 



M.K(o-(0))/mp G V(u ,e) and 
m Ajv (cr(t)) ^ - 3e 



(3.24) 



(3.25) 
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where the second inequality is a consequence of Markov's inequality. Now we fix k G N* and 
a > and call 



C 



M K ( ({i + l)t 
a 



mp 



M 



K 



mp 



a 



it 



^ e, yi < k 



L 1 



the event of continuity, and 



Mi ^ k, 3vi G BV a : 



M 



K 



mp 



a 



it 



G V{vi,e) 



the event that the magnetization profile is close to BV a at any time it/k. Then we remark that 
£K(a(0))€V(«o,e), 



m AN (a(t)) ^ mp - 3e 
and C n L> a occurs 



3v = (vi) i=0 _ k e c 8e/m/3 (u ), 
^Kf))ev(«,4 



(3.26) 



Indeed, the definition of D a yields a sequence t> = (fj)j=o...fc- We can take t>o = since we know 
that M.k ip (0)) / m /3 ^ V(no,e). This sequence is 3e-continuous according to the properties of 
the Vi and to the definition of C. Finally, we have 



I 1 - wjfcllii 



[0,1]° 



(l-v k (x))dC d (x) 



(1 - M K (cr(t))(x)/mp) dC d (x) + e 

[o,i] d 

^ l-mn N (a t )/mp + 2e 
^ he/mp 



(3.27) 



and therefore Mk (p (t)) /mp G V(l, Qe/mp). So we can fix v k = 1, which makes the evolution 
8e/m^-continuous. 

Now we conclude the proof of (3.24). As a consequence of (3.25) and (3.26) we have the 
inequality 



ep J ^ P 



J,A N ,+ / 3-U = (wj)i = o...A: G C S£ / mi3 (u ), 



Pa* 



and therefore, according to the invariance of f° r the Glauber dynamics, and to the definition 
(3.17) of the bottleneck set £> e (iio,r r ), 



P J ^ p{+pi+ pl 



where 



Pi 



Pi 



4 



A" 



m/3 
m/3 



gH 



G* V(BV a ,e) 



Ma- <t 



<<iv 



-Ma-KO)) 



L 1 
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Finally we bound each contribution separately. First, according to Lemmas 3.15 and 3.17 
have 



we 



limP ( a J : + 



M 



K 



B 8e ) < e -N d ~\lC r {u )+^{u )-20 



i=l. ..n. 



1 



for e > small enough, for some 7 > 0, and for any (uq, t t , 6, 7)-covering (TZi)i=i... n of 8*uq 
with 5 > small enough, for any K large enough. So for 

k = kN = [exp(N d ~ 1 (K7(uq) — £/2))] 

we have, under the same conditions, 

1 



iv 



6 



(3.28) 



Then, the exponential tightness property (see (3.18) above) implies that, for a > large enough, 



limP ( P2 ^ -e 



v 



^ d - 1 (J"'-( U o)-€) 



6 



j / i=l. ..n . 



and finally Lemma 3.19 implies that, for large N, uniformly over J and t ^ -Kjv> 



pi < exp -ciV a . 



Summing the last three displays yields (3.24). 



(3.29) 

(3.30) 

□ 



3.3.5 Proofs of Theorems B31 and 12151 



Here we conclude the proof of Theorems 2.2 



and 



averaged autocorrelation A^(t) defined at (1.27) to the dynamics of the overall magnetization 
m\ N defined at (3.16). 



2.5 We state one more Lemma that relates the 



Lemma 3.20. For any e > 0, N € N* and A ^ 1 one has 



A A {t) > e 2A E 



(3.31) 



Proof. Minkowski's inequality implies, for any v ^ 1, that 

W J \T J (t)n 9 -n J '+ir x \ v dn 



|T J (t)m Ajv - ^' + {m kN )\ v 



l/u 



l/u 



where ir x : £ — » P is the function which associates, to the spin configuration a G S, the spin at 
x, <r(x). Taking f = 1/A, the translation invariance of E and of the Glauber dynamics implies 
that 



A x (t) > E 



\T J {t)m hN - /'+(m Ajv )| 2 



29 



Hence, for any e > we have 

A x (t) > £ 2A E [l W , +{mAN)>m ^ £} x (//'+ {T J (t)m AN <mp- 2s)) 
and we conclude by using the attractivity of the Glauber dynamics: 

//>+ (T J (t)m AN ^ mp - 2s) ^ /'+ (r/^(t)m Ajv < - 2e) 



(3.32) 

□ 



Proof. (Theorem 2.2). Let (5 > and A ^ 1. We assume that <Y\ < oo, otherwise there is 
nothing to prove. We fix (uq,t v ) G IC and £ G (0,/C r (ito)) such that 



1 + < * A + <ty2. 



Then, for any t > we call iV(t) the smallest integer A" such that t ^ exp(A^ rf_1 (JC r (uo) — £))• 
According to Lemma 3.20[ to Propositions 3.14| and |3.18| > to Lemma 3.10[ for any s > small 
enough we can find 6, 7 > such that, provided that N(t) is large enough, 

A x (t) > ^exp^A^-^J^M + O-^W^H^M + O)- 
The definition of N(t) implies finally that 

' , A (T r (u ) + e) + (I r («o) + / iV(t) x 



A A (t) > -exp 



-(logt)- 



AC(uo)-C 

for i large enough, and this gives the claim as A r (t) — > +00. 



iV(t) - 1 



□ 



Proof. (Theorem 2.5) Let 5 > and assume that k > 0. There exists (uo,T r ) G IC and £ > 
such that 

Given A" G N* large, we define an intermediate side-length M = M(N) as the largest integer 
M such that 

N d > exp(M d - 1 (X r (n ) + 2e)). 

First, we prove that the event of dilution, on scale M, occurs with large probability inside Am- 
More precisely, we consider L = M + [-i/M] and 

X N = {0, . . . , [N/L] - l} d . 

We also consider 5 > and a 5-covering for duo, and for any x G AT we call 



<M 



{dilution occurs in Lx + Am} 
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the Lx-translate of Qm- According to Lemma 



3.10 



we have, for 5 small enough, 



p ( n (^r ) 

\x&X / 



< (l-exp(-M d - 1 (X'-K) + 0) 



[N/L}° 



cxp 



N 

T 



exp(-M d - 1 (r( Uo ) + 6) 



^ exp 



1 



(7 log A 



exp(M^ 1 



^ exp -N 2 ( xr ( u o)+ 2 



(3.33) 



where we use the definition of M and the fact that L ^ Clog A, for some C < oo. We call 
now Vf^ = P| C\ zeX z <x (^m) C ^ ne even t that x is the smallest box for which dilution occurs, 



where X is ordered according to the lexicographic order. According to Proposition 3.18 and to 
the attractivity of the dynamics, the conditional probability 



P {An { T An (^(M^ 1 (/C> ) - 6)) ™ Lx+Am < m p - 2s) > exp {-M d ~ l (T r (u ) + 0) 



V 



goes to 1 as A — > oo (and thus M). On the other hand, Proposition 3.14 and the monotonicity 
of the system imply that, as well, the conditional probability 



P ( Alt ( m Lx+A M ) ^mp-e) 



T) x 



M 



goes to 1 as A — > oo. Therefore, we have shown that the P-probability that there exists x £ Xjy 
with both 

A + N (?a£ ( ex P( Md_1 (£>o) - 0)) ™ Lx+Am ^ mp - 2e) ^ exp {—M d ~ 1 (T r (u ) + £) 



At ( m -C*+Ajvr) > rnp 



e 



goes to one as N — > oo. Taking tulx+Am as a test function in ( 1.30 ) proves the lower bound on the 
relaxation time. For obtaining the lower bound on the mixing time one can consider / = mi jX+AM 



in (1.32), and oq any initial configuration for which T A ' + (exp(M ci 1 (KZ r (uo) — £))) m Lx+Am ^ 



mp — 2e, which is possible as this event has positive probability. 



□ 



3.4 The geometry of relaxation 



The first part of this Section is dedicated to the proof of Proposition 2.1 Then we show that the 
gap in surface energy can be computed on more restrictive evolutions, and finally we compute 
IC r (uo) for two simple initial configurations. 



3.4.1 A lower bound on the gap in surface energy 



The statement of theorems 2.2 and 2.5 is non-empty only if K7{uq) happens to be strictly positive 
for some initial profiles {uq,t t ) £ IC We state the next Theorem, which has Proposition 2.1 as 
a Corollary: 
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Theorem 3.21. Let (uo,T r ) E IC. Assume that the boundary of uq is C . There exists a non- 
decreasing function H : R + — > R + ; with H{5) > 0, \/5 > 0, that depends only on (uq,t t ), such 
that 



JC r {u ) ^ sup r q {n) x H 



max(sup ne5 d-i Ti{n), 1) 



(3.34) 



Proof. (Proposition 2.1). When r mm (n) < r 9 (n) for all n E S and the boundary of no is 
C , it follows at once from (3.34) that IC r (uo) > 0. Now we address the proof of (2.2). We 
discuss first the consequences of the assumption P( J e = 0) > 0. When J = along a section of 
a rectangle 1Z the surface tension in that rectangle is = 0, thus for all e > 0, all /3 > 0, 



I n (e) < -||n||ilogP(J e = 0). 
Now we consider the initial configuration defined by 

= Xb(z ,i/4) and t t {x) = min(l, r <? «°)/2), Vx E d*u 



(3.35) 



where zo = (1/2, . . . , 1/2). According to (3.35) there is C < +oo not depending on j3 such that 
the cost of dilution I r (uo) defined at (1.22) satisfies X r (uo) ^ C. We also have F t {uq) ^ C 
where J- v (uq) is the initial surface energy defined at ( |1.20 ), and C is the perimeter of uq. So we 
conclude already that the numerator X t {uq) + AJ-" r (tto) in the definition (2.1) of X\ is bounded 
by a constant which does not depend on (3. Finally we recall that, thanks to the assumption 
that P(J e = 0) < 1 - p c (d) and to Proposition 2.13 in |Wou09j . 



lim inf inf 

/3->oo ne5 d - 1 



rj{n) 

P 



> 0, 



(3.36) 



and therefore the inequality (3.34) implies that, for some c > and for all fH large enough, 



fC r (u ) > 



1 inf r 



6S d-ir«(nj 



sup T q (n) x H - 
ieS d-i \2swp neS d-iTi(n) 



The convexity and the lattice symmetries of the Wulff crystal W q imply that the ratio 

is bounded from below by some constant Cd not depending on /3, hence 

K, r \u ) > sup T q (n) x H{c d ) ^ CP 
ries 11 - 1 



for large j3. The claim (2.2) follows. 



□ 



proof. We introduce the new surface energy 



The proof of Theorem |3.21| requires three Lemmas, that are stated now together with their 

T r dH d ~ l , VueBV. (3.37) 



d*u\d*UQ 



T q {n u )dU 



d-l 



d*und*uo 
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Lemma 3.22. Let (no = Xu , tT ) £ IC and u = xu £ BV. Then, the profiles v = Xu \u an d 
w = Xu\u satisfy 



Proof. We remark first that 



F r >-{v)+F r >-{w). 



f T q {n u )dU d - 1 - f 



T r dH 



d-l 



d*uo\d*u 



(3.38) 



(3.39) 



As stated in Theorem 3.61 in [AFPOOj . given any u = xu G BV the local density of U at x is 
either 0, 1/2 or 1, for "W^-almost all x G R , and the set of points at which the local density 
is 1/2 coincides, up to a % d_1 -negligible set, to the reduced boundary d*u (denoted FU in 
[AFP00]). This implies the two equalities 



d*u\d*u = {d*v \ d*u ) U (d*w \ d*u ) , 
d*u \d*u = {d*wr\d*u )u(d*vr\d*u ) 



(3.40) 
(3.41) 



up to H d ~ x -negligible sets, where U stands for the disjoint union (again, up to % rf-1 -negligible 
sets). Furthermore, the outer normal at x G d*w \ d*uo (resp. x G d*v \ 8*uq) corresponds 
(^ o!_1 -a.s.) to the outer normal at x G d*u\d*uo (resp. to the opposite of the former). Thus, in 
conjunction with (3.39), equations (3.40) and (3.41) imply respectively the r 9 and the r r part 
of ( pT38l >. □ 

Lemma 3.23. Let (uq,t v ) G IC and u = xu 6 BV. For any h > 0, there is a finite collection 
of droplets (i>i = XVi)ieI such that 

i. The Vi are disjoint and their union is U , 

ii. Each Vi is included in a box of side-length at most h, 
Hi. And 



F r 



u 



iel 



h 



(3.42) 



Proof. We define I = {0, . . . , \l/h] } and call 

Vi = {x G U : ^ Xk — Zk — ikh ^ h, Vfc = 1, . . . , d}, \/i G I 



where z is some point of (0, h) d . The first two properties hold trivially. Now we chose z so that 

^({xGU :(x k - Zk )/h€N})dz k = C d (U) 



(3.42 ) holds. For every k G {1, . . . , d} we have 

f 

Jo 



and therefore we can chose z k G (0, h) such that H ({x G U : {x k — z k )/h G N}) ^ £ d (U) /h = 
||1 — u\\li /(2h). Equation (3.42) follows as the new portions of interface in the Vi that were not 
present in v are exactly the {x G U : {x k — z k )/h G N} as k = 1, . . . , d. □ 
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Lemma 3.24. Assume that u G BV, no G S d 1 and e > satisfy 

T r (x) + {e + \\n u x -nv\\ 2 ) sup r«(») < T q (n ) (3.43) 

/or a// xed*un d*u \ N, where H^iN) = 0. Then 

F r -(u) > Cd £SupT«([|l-u|| L i)^ (3.44) 
where Cd > is a constant that depends only on d. 
Proof. We let 

W q = {zelR d : z-n^T q (n),\/ne S 11 - 1 } 
and W = {z G W 9 : z • n£ -r r (x),Vx G d k und k u \M^. 

The surface tension associated to W is 

f(n) = sup z • n, Vn G S 1 ^ 1 

and it satisfies obviously the relations 

f(n) < T 9 (n),Vn G 5 d_1 
f«) < — T r (x), G <9*« n <9*iio \ A/" 

therefore 

r>-(u) ^ J f{n u x )m d -\x). (3.45) 
The isoperimetric inequality for f (see |KP94j ) implies in turn that 

r''-(u) > dC d {W) 1/d ( ^-L^lA " (3.46) 



so the claim will follows from an appropriate lower bound on the volume W. There exists 
zq G W q such that zq ■ uq = r 9 (no). Now we take z G W q , close to — zq, and prove that it lies 
in W: if x G d*u n d*u \ M, then 

z ■ ri% = (z + zq) • n" - z ■ (n" - n ) - z • n 

< Ik + zoh + sup r 9 ||n" - n ||2 - T q (n ) 

S d-i 

as the Euclidean norm of zq is not larger than r^zo/H^olh) ^ sup S d-i r 9 . According to our 
assumption (3.43) it follows that 

W D W q D B[-z ,e sup r q 

The set W 9 is convex and contains all the images of zq by the symmetries of therefore it 
contains an hypercube which vertices are the images of zq by the above mentioned symmetries. 
Consequently the volume of VV is at least l/2 d of the volume of the ball B (— zo^supr 9 ) and 
the proof is over. □ 
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Proof. (Theorem |3.21 ). Given 5 > 0, the continuity of r r and the smoothness of du$ imply that 
there exists h{5) > such that 

As we claimed, h{5) depends only on (uo,T r ), and of course on 5. Also, it is clear that one can 
chose the function h non-decreasing. Furthermore, to a given direction n G S^ 1 corresponds 
z G W q such that z ■ n = T q (n) and therefore 



T q (n') > z-ri 
> r q (n) 



|n — n'W? sup r q . 



It follows that any u G BV with diameter not greater than h(5) satisfies assumption (3.43) of 
Lemma 13.241 when one takes 

£ = 5 mU d , UQ (r q (n^)-r r (x)) _ 
6 max(sup ng5 d-i r 9 (n), 1) 



Now we consider some phase profile u G BV, call v and w the profiles given by Lemma 3.22 and 
then (vi)i£j the union of the droplets associated to both v and w by Lemma 3.23 In conjunction 
with Lemma 13.241 we obtain 



F( U )-r\u Q ) > C^SUpT^CIl-^llix)^ 

iei 



1 — v\\ L i + ||1 — w\ 



Li, 



h(S) 



According to the definition of v , w and of the Vi we have 

II^-^IIl 1 = 



II 1 - v \\Li + II 1 - W\\ L 1 



therefore, the trivial inequality x a + y a ^ (x + y) a for a G (0, 1) and x,y > implies that, for 
every u G BV, 



.^(u) - ^(no) ^ c rf (5supr 9 (||n - u ||li)" 



dT q ( ei 



But x 1 1 / a! — bx is maximized for . 



< d-l s 
\ bd i 



and equals c d (l/b) d 1 . So, taking b = £^^sK(S)' 



for any u G BV at L 1 -distance (^r) from u (this distance is smaller than the L 1 distance 
between uq and 1 if we require that sup,5 >0 h(5) be small), we get 



F(u) -F{uq) > c d 5 sup r q {5h(5)) d - 1 
and we have proved (3.34) with H(65) = d d 5 d h d ~ 1 (5), as e = 5. 



□ 
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3.4.2 Continuous evolution and continuous separation 

In this paragraph we give an alternative formulation of the surface energy gap K7{uq). We 
introduce a set of continuous evolution of the phase profile, for which the boundary splits 
continuously from its original location: 



C(u ) 



v : t G [0, 1] i-)- Vt G BV: vo = uo, v\ = 1, 
t 1 — y Vt is continuous for the L 1 -norm 
and t 1 y 1-d*v t r\d*u is continuous for the 
L 1 -norm associated to the measure 1~L d ~ 1 



(3.47) 



and then define /C£ s , the gap in the free energy associated to the optimal droplet removal in this 
class of evolutions: 

V cs (u )= inf sup F r (y t ) -J^{uq). (3.48) 
oeC(wo) te[o,i] 

Theorem 3.25. For all {u ,T r ) G IC, 

K, r cs (u ) = K, r (u ). 

The inequality /C£ s (ito) ^ lC r (uo) is clear since the set C(uq) is more restrictive that those 
evolutions considered in the definition (1.24) of K7(uq). We prove the reverse inequality with the 
help of an interpolation: given an evolution v = (i>i)j=o...fc £ C £ (uq) (see (1.23)) we interpolate 
from v a continuous evolution v' with continuous separation from d*uo, at the price of a small 
increase in the maximal cost, negligible as e — > (and uniform over v): 

V G C(uq): v s iu = Vi and sup J^iv't) ^ max J^ivi) + o (1). 
1 te[o,i] i=0 - k 



The next lemma is one of the keys to the proof of Theorem 3.25 



Lemma 3.26. Let (uq, t t ) G IC and 5 > 0. For any Borel set A C [0, l] d with volume £ d (A) ^ 
8, there exists a collection of measurable sets U = (Ut)tem. such that: 

i. 1 1—)- Ut is a non- decreasing function with 

lim Ut = $ and lim Ut = T 

t— >— CO t— >+oo 

where T is the tube T = [0, l] d x R. 
ii. The function 1 1-> Ut — te^ is 1-periodic. 

Hi. The volume 1 1— > C d {Ut fl [0, l] 01 ) zs a continuous function of t 

iv. The area uo) is a continuous function oft 

v. The portion of the boundary of Ut that intersects A + Ze^ in T has a small area: 



sup U d ~ l (dUt n (A + 7Le d ) nf K7\4 

te[o,i] v y 



for 5 > small enough. 
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Proof. See Figure [2] for an illustration of the proof. To begin with, we partition [0, l] d in 
horizontal slabs: let n = \ \/y/5\ and call 

,j i — 1 i 

x G 0, If : x ■ e d sC - 

n n 

for i G {l,...,n}. Because each slab has a volume at least V~5, for each i G {0, . . . n — 1} there 
exists Zi G (z/re, (z + l)/n) such that the density C d ~ l (A n {z : z ■ e d = Zi}) of A at height Zi is 
not larger than \J~5. We extend then the definition of zi by periodicity, letting 

z i+n = 1 + Zi, Vi G Z. 

Then, we let n = cosaei + sinae^ for some a G [0,7r/3] such that d*uo has no face orthogonal 
to n and define 

U t = \xGT: ^fW" j VteR 

[ a; • e d ^ z\ nt -\ and x ■ n ^ k J 

where / is the piecewise linear function defined by: Vi G %, 

l{i/n)+ = z i e d ■ n 

l {i+i)/n = (zi+ie d + ei) • n 

and I linear on each interval (i/n,(i + l)/n]. The set Ut evolves as follows: between times i/n 
and (i + l)/n, Ut invades the region {x G [0, l] d : z% ^ x ■ ^ ^i+i} by the mean of a front line 
normal to n, that moves at a constant speed. 

It is immediate from the definition that Ut — ted is 1-periodic and that Ut increases contin- 
uously in volume. The 'H d ~ x measure of Ut H d* uq is non-decreasing and the assumption on n 
ensures that it increases continuously. We consider at last the portion of the surface of Ut in 
T that might intersect A + 7Le<i- We just have to take into account the upper portion of dUt, 
made of the two planes at height z%, -Zj+i, and of a portion of plane normal to n. Recall that the 
Zi have been chosen so that the density of A + Ze^ at height z% does not exceed VS. Similarly, 
because a ^ tt/3, the piece of plane orthogonal to n has a surface at most 4/n ^ 5/V5 for 5 > 
small enough. The claim follows. □ 

The second key argument is periodicity: as seen on Figure [2j the interpolation between vo 
and v\ has to choose first the region where the cost of v\ is smaller than that of v$ — which 
means that we have to fix to carefully. 



Proof. (Theorem 3.25). Let 5 > 0. There exists e G (0, 25] and fcgN, together with v G C £ (uq), 
such that 

max J^M < KJiua) + 5. 

i=0...k 

Starting from v, we construct a continuous evolution u G C(uq) that has a maximal cost not 
much larger than that of v. It is enough to do the interpolation between two successive Vi, as 
one can paste together the successive interpolations to deduce the continuous evolution u. 
Hence we consider vq,v\ G BV and assume that \\vq — ^ 25. We let 

A = {x : vq(x) / vi(x)} 
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Figure 2: The construction of Ut, and the interpolation ut between vq and v\. 



which has a volume at most 5. Lemma 3.26 applies and there is Ut with properties (i)-(v). Given 
to G R and t € [0, 1] we let 

G t(h t = {x£ [0, l] d :3keZ: x + ke d £ U to+t \ U to } , 

for any to the set Gt ,t increases continuously from the empty to the full set in [0, l] d , makes the 
surface T~L d ~ 1 (d*uo fl G^t) a continuous function of t, and the area 

u d - x (dG t0 ,t n A n r) < 14v^ 

small, for 5 > small enough. Now we define 

„ /_i _ / M x ) if z i G to ,t 
u ^ x >-\ Vl ( x ) ifxeG t0 , t . 

The cost of ut decomposes in the following way: it is the sum of the cost of vo in G£ t , of the 
cost of u\ in Gt 0t t, and of the cost of dGt Q t in A. In other words, 



r\ Ut ) ^ r\vo) - r G (vo) + r Gtn » + 



(3.49) 



where ^^(u) stands for 



T r dn d - 1 



nd*u r\E 



-r / T q (n u )dH 

{d*u\d*u Q )r\E 



d-l 



It is clear that the initial cost of Ut is J- r (vo) and that its final cost is J- r {vi). Yet in the interval 
(0, 1) it could be that Gt 0: t selects first the region where v\ has a larger cost than vq, leading 
to a maximal cost larger than expected. We rule out this possibility with an appropriate choice 
for to - see Figure [2] for an illustration of the discussion below. 
For to £ R and t £ [0, 1) we consider 



f{^t) = r Gtot ( Vl )-r GtQt {vo). 
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Our aim is to extend / to arbitrary values of t G E + . For k G Z, we denote by the translated 
of by ke^, then for any to G H and t E [0, 1) we have 



fees 

from the definition of Gt ,t- The latter formula permits to extend / to R x R + and puts in 
evidence the existence of a function g : 1R — >• 1R such that 

/(to, *) = ff(to + 1) - <?(to), V(t , t) G E x R+. 
This function is, apart from a linear correction, 1-periodic: for all t G E, 

g(t+l)=g(t)+T r (v 1 )-T r (v ), 

in other words, 

g{t) = h{t)+t{r(v l )-F{v Q )) 
where h is a 1-periodic function. Now we fix to such that 

h(to) ^ sup h(t) — \/5, 
t 

it is immediate that 

f(t ,t) = g(t + t)-g(t ) 

= h(f +t)-h(t )+t(r r (v 1 )-r r (v )) 
^ *(.t> 1 )-.f>o)) + y / £. 



Reporting into (3.49) we conclude that ut is a satisfactory interpolation between vq and v±: 



provided that 5 > is small enough, 

F{u t ) < (1 - t)^ (w ) + tT (t>i) + 15v^, Vt G [0, 1). 

□ 

3.4.3 The gap in surface energy in the isotropic case 



In this paragraph we use Theorem 3.25 to compute the gap in surface energy in the two dimen- 
sional, isotropic case. 

Lemma 3.27. Assume d = 2, T q (n) = 1 for all n G S^ 1 and consider 

u o = Xb and r r (x) = A, Vx G d*uo 
where B is the disk of radius r < 1/2 centered at (1/2, 1/2), and A G (0, 1). Then 

JC r {uo) = 2r \/l — A 2 — A acos A . 
An optimal continuous evolution in this setting is illustrated on Figure [3j 
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Figure 3: A continuous evolution of minimal cost. 



Proof. In order to simplify the notation we will consider r = 1/2. The upper bound 

K7(uq) ^ sup {sm.6 — \9) 

6»e[0,7r] 

is immediate if one considers the continuous evolution (ut)te[0,l] defined by 

, . f — 1 if x G B and x ■ e2 ^ 1 — t 
Ut{x) = { 1 else 

and 6 satisfying 1 — t = 1/2 + sin9, as illustrated on Figure [3| The lower bound is scarcely 
more difficult to establish: given (ut)$e[o,i] £ £(^o) a continuous evolution with continuous 
detachment, there is t £ (0, 1) such that 

^(d^uo \ d*u t ) = acos A. 




Figure 4: Reduction of u to a portion of disk in three steps. The surface energy decreases, the 
length of contact is preserved. 

Optimizing the droplets of ut as in Figure [4] - we replace each portion of the interface not 
in d*uo with a segment (step (1)) - we obtain a profile u' t with a lower cost, yet it still has the 
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same contact length acos A with 8*uq . By isotropy of surface tension it is possible to aggregate 
the droplets together (step (2)) and obtain u'( with a unique droplet and a lower cost, preserving 
again the length of contact. At last, inverting the order of the segments and arcs and optimizing 
again we see that the profile of lower cost that satisfies TV-^uq \ d*u) = acos A coincides, apart 
from a rotation, with the profiles considered in the upper bound. The claim follows. □ 

3.4.4 The gap in surface energy when the Wulff crystal is a square 

Here we compute the gap in surface energy for another simple case, when the Wulff crystal W q 
is a square. We also show that the gap in surface energy is strictly bigger than the cost of the 
less likely overall magnetization. 

Lemma 3.28. Assume d = 2, T q (n) = \\n\\i for all n G S^ -1 and 

u o = Xc an d t t (x) = A, Vx G d*uo 
where C = [1/2 - r, 1/2 + r} 2 , r < 1/2 and A G [0, 1]. Then, we have 

K7{uq) = 2r[l-A]. 

Proof. Again we consider r = 1/2 in order to simplify the notations. The upper bound on the 
additional cost is immediate considering ut = X{x-e 2 ^i-t}- F° r the lower bound we need a finer 
analysis. First, it is a consequence of the assumption on r q that for any open, connected PcE 2 
with extension h x ,h 2 in the canonical directions, that is: 

h(k) = supx • — inf x • e^, 

we have 

F q ( X u) >2h 1 + 2h 2 . 

Then, we decompose a profile configuration u into its droplets (C/j)j^o- We call hj,hf the 
extension of U\ in the canonical directions and let ij the length of contact between d*xiii an d 
9*uo, so that, for all i: 

r-( XUi ) > 2/i 1 + 2h 2 - (1 - X)h 

If a droplet Ui touches two opposite faces of d*uo, say h 1 = 1, then its extension in the orthogonal 
direction is at least h? ^ (U — l)/2 and the inequality 

follows. If on the opposite the droplet is in contact with at most two adjacent sides of d*uo, we 
have h 1 + h 2 ^ k and hence 

^(Xu,) > (l + X)k. 

Assume now that the total length of contact is 3, i.e. that J2i k = 3. A consequence of the 
former lower bounds is that, whether or not some droplet touches two opposite faces, the cost 
of u is at least ^(u) ^ 1 + 3A. The claim follows. □ 
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We conclude this work with a comparison between the bottleneck due to the positivity of 
JC r (uo) and the one due to the continuous evolution of the magnetization: 

Lemma 3.29. In the settings of the former lemma, with furthermore A 6 (0, 1), we have 

K7{uq)> sup inf F r {u) — F 7 " {uq) . 

me[-l,l] ueBV: J[o,i]d u=m 




Figure 5: The two profiles u x m and v? m : for m ~ 1 the first one is better, for m ~ — 1 the second 
one has a smaller cost. 

Proof. We provide an upper bound for the right hand term, considering for a given m the two 
profiles (see Figure [5]) 

x _ / -1 if x e [0, l] d and x • e 2 < ^ 
1 else 



and 



u 



2 _ J -1 if x G [0, l] d and min(x • ei,x- e 2 ) < 1 - 



[ 1 else 

that both satisfy the volume constraint f, Q ^ u = m. It is immediate that 



F 7 \u l rn ) = 1 + (2 - m)X and F(u 2 m ) = 4A + 2(1 - A) ' ' " "' 
hence 

sup inf F r '(u) ^ sup min (jF(u m ),'F(u% n )') . 

me[-l,l] -"GBV:/ [0|1]d «=m me[-l,l] 

Note that F{v}^) decreases with m while .F(ii^J increases with m. Because of their extremal val- 
ues there exists some mo £ (0, 1) at which /'(ii^J = /(u^J = sup me r j^i min (^(it^J, ^(u^)) , 
and since mo < 1 we have in particular ^(ti^J < J r (uL 1 ) = 1 + 3 A which is the maximal cost 
of an optimal continuous detachment evolution. □ 
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